J. Pure Appl. & Ind. Phys. Vol.2 (1), 55-62 (2012) 



The Onset of Convection in Walters' (Model B ) fluid in a 
Darcy-Brinkman Porous Medium 

G. C. RANA 

Department of Mathematics 
NSCBM Govt. P. G. College, Hamirpur-177 005, Himachal Pradesh, India 

ABSTRACT 

The onset of convection in Walters' (Model B ) elastico-viscous 
fluid in porous medium heated from below is considered. For the 
porous medium, the Brinkman model is employed. By applying 
normal mode analysis method, the dispersion relation has been 
derived and solved analytically. It is observed that the medium 
permeability, gravity field and viscoelasticity introduce oscillatory 
modes. For stationary convection, the Darcy number has 
stabilizing effect whereas the medium permeability has 
destabilizing effect on the system. The effects of medium 
permeability and Darcy number has also been shown graphically. 

Keywords: Walters' (Model B) fluid, thermal convection, 
viscoelasticity, porous medium. 



1. INTRODUCTION 

In recent years, considerable interest 
has been evinced in the study of thermal 
instability in a porous medium, because it 
has various applications in geophysics, food 
processing and nuclear reactors. A detailed 
account of the thermal instability of a 
Newtonian fluid, under varying assumptions 
of hydrodynamics and hydromagnetics has 
been given by Chandrasekhar 1 . Lap wood 4 
has studied the convective flow in a porous 
medium using linearized stability theory. 
The Rayleigh instability of a thermal 
boundary layer in flow through a porous 
medium has been considered by Wooding 11 . 



There are many elastico-viscous 
fluids that cannot be characterized by 
Maxwell's constitutive relations or 
Oldroyd's constitutive relations. One such 
class of fluids is Walters' (Model B ) 
elastico-viscous fluid having relevance and 
in chemical technology and industry. 
Walters 10 reported that the mixture of 
polymethyl methacrylate and pyridine at 
25°C containg 30.5g of polymer per litre 
with density 0.98g per litre behaves very 
nearly as the Walters' (Model B ) elastico- 
viscous fluid. Walters' (Model B ) elastico- 
viscous fluid form the basis for the 
manufacture of many important polymers 
and useful products. A good account of 
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convection problems in a porous medium is 
Vafai and Hadim 9 , Ingham and 

.5 



by 



given 

Pop 2 and Nield and Bejan 3 



The investigation in porous media 
has been started with the simple Darcy 
model and gradually was extended to Darcy- 
Brinkman model. Kuznetsov and Nield 3 
have studied thermal instability in a porous 
medium layer saturated by a nanofluid: 
Brinkman model. Sharma and Rana have 
studied thermal instability of a Walters' 
(Model B ) elastico-viscous in the presence 
of variable gravity field and rotation in 
porous medium. Recently, Rana and Kango 6 
have been studied the effect of rotation on 
thermal instability of Compressible Walters' 
(Model B ) elastico-viscous fluid in porous 
medium Rana and Kango 7 have also studied 
thermal instability of compressible Walters' 
(Model B ) elastico-viscous rotating fluid 



permitted with suspended dust particles in 
porous medium 

The interest for investigations of 
non-Newtonian fluids is also motivated by a 
wide range of engineering applications 
which include ground pollutions by 
chemicals which are non-Newtonian like 
lubricants and polymers and in the treatment 
of sewage sludge in drying beds. Recently, 
polymers are used in agriculture, 
communications appliances and in bio 
medical applications. Examples of these 
applications are filtration processes, packed 
bed reactors, insulation system, ceramic 
processing, enhanced oil recovery, 
chromatography etc. 

Keeping in mind the importance in 
various applications mentioned above, our 
interest, in the present paper is to study the 
onset of convection in Walters' (Model B') 
fluid in a Brinkman porous Medium. 



g = g(0,0,-g) 




Porous Medium 



rTTTTTT 



Walters' (Model B) 
fluid 



X 



Heated from Below 



Schematic Sketch of Physical Situation 
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2. MATHEMATICAL MODEL AND 
PERTURBATION EQUATIONS 

Here we consider an infinite, 
horizontal, incompressible Walters' 
(Model B ) elastico-viscous fluid of depth d, 
bounded by the planes z = 0 and z = d in an 
isotropic and homogeneous medium of 
porosity s and permeability k 1; which is 
acted upon by gravity g(0, 0, -g). This layer 
is heated from below such that a steady 

adverse temperature gradient B |^|) is 

maintained. The character of equilibrium of 
this initial static state is determined by 
supposing that the system is slightly 
disturbed and then following its further 
evolution. 

Let p , v , v ' , p, G, T, p., a and v(0, 0, 0), 
denote, respectively, the density, kinematic 
viscosity, kinematic viscoelasticity, pressure, 
medium porosity, temperature, effective 
viscosity, thermal coefficient of expansion 
and velocity of the fluid. 

The equations expressing the 
conservation of momentum, mass, 
temperature and equation of state for 
Walters' (Model B ) elastico-viscous fluid 
are 

Mf + ;c-H:-> + «( i+ S)- 

±(v-v'^-)v + ^V 2 v, (1) 

fci V dtj p 0 

V.v = 0, (2) 
E^+(v.V)T = /cV 2 T, (3) 

p = p 0 [l-a(T-T Q )], (4) 

where the suffix zero refers to values at the 
reference level z = 0. 
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Here 

E = £ + (1 - £) (^) 

\PoCfJ 

which is constant, k is the thermal 
diffusivity, p s c s ; p 0 cy denote the density 
and heat capacity of solid (porous) matrix 
and fluid, respectively. 
The initial state of the system is taken to be 
quiescent layer (no settling) with a uniform 
particle distribution number. The initial state 
is 

v = (0,0,0), T = -Bz + T 0 , 

p = p 0 (l + aBz), (5) 

is an exact solution to the governing 
equations. 

Let v(u,v,w), VdQ.,v,s), 9, Sp and Sp 
denote, respectively, the perturbations in 
fluid velocity v(0,0,0), the perturbation in 
particle velocity v d (0,0,0), temperature T, 
pressure p and density p. 

The change in density Sp caused by 
perturbation 8 temperature is given by 

Sp = -ap 0 Q. (6) 

The linearized perturbation equations 
governing the motion of fluids are 

1 dv 1 _ r Sp 1 
= VSp — q — 

E dt Po Po fei 



(v-v'-pjv + ^V 2 v, (7) 

V dt) p 0 

V.v = 0, (8) 

(7TL d . . \ 

--+l)v d =v, (9) 

E d ± = Bw+ kV 2 8, (10) 



where w is the vertical fluid velocity. 
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(V*w) = -i(„-„'£)v*w 



In the Cartesian form, equations (7)-(10) 
with the help of equation (6) can be 
expressed as 



+ 



(17) 



ldu Id 1 



u + -V 2 u, 

Po 



(11) 



Idv Id 1 
a ot po oy k-i 
(v — v'^-) v + — V 2 v, 

V dtj po 

ldw 1 d 1 

-— = — (Sp) - gaB - — 

e ot po oz /q 

(v - v'^-) w + — V 2 w , 

v at) p 0 



du dv dw _ g 

5% dy dz 



(12) 



(13) 
(14) 
(15) 



Operating equation (11) and (12) by — and 
^ respectively, adding and using equation 
(14), we get 



i a fdw\ i f- 2 d 2 \ _ 



where <" = ^ - ^ is the z-component of 
vorticity. 



where V 2 = £_ + ^_ + £_ 

dx 2 dy 2 dz 2 



3. THE DISPERSION RELATION 

Following the normal mode analyses, we 
assume that the perturbation quantities have 
x, y and t dependence of the form 

[w,s,0,(] = [W{z),S{z),6{z),Z(z)\exp(ilx + 
imy + nt), (18) 

where / and m are the wave numbers in the x 

1 / 

and y directions, k = (I 2 + m 2 ) '2 is the 
resultant wave number and n is the 
frequency of the harmonic disturbance, 
which is, in general, a complex constant. 
Using expression (18) in equations (17) and 
(15) become 



d 2 
dz 2 

2 



W = -gk 2 a8 - —(y-v'ri) 
k 1 



EnQ = pW + k - fc 2 ) Q. 



(20) 



Equation (19) and (20) in non dimensional 
form, become 



[| + ^] (D 2_ a 2 )M/ + ^!^ + 



Operating equation (13) and (16) by £ (D 2 - a 2 ) 2 W = 0, 

(v 2 — and respectively and adding 

to eliminate 8p between equations (13) and 
(16), we get 



[(D 2 - a 2 ) - E Pl a]8 = - P —W, 



(21) 
(22) 



where we have put 
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a = kd, a = — , F = — and P ; = — , is 

v d 1 1 d l 

the dimensionless medium permeability, 
p x = - , is the thermal Prandtl number. 



Eliminating 0 between equations (22) and 
(21), we obtain 

[l + (7 - F) <t] (D 2 - a 2 )(D 2 - a 2 - 
Ep l( r)W - Ra 2 P t W + D A (D 2 - a 2 ) 2 (D 2 - 
a 2 -E x y x o~)W = 0, (23) 

where R = 3<x ^ d [ s the thermal Rayleigh 

number and D A = ^% , is the Brinkman- 

A fid 2 

Darcy number modified by the viscosity 
ratio. 

Here we assume that the 
temperature at the boundaries is kept fixed, 
the fluid layer is confined between two 
boundaries and adjoining medium is 
electrically non-conducting. The boundary 
conditions appropriate to the problem are 
(Chandrasekhar 1 ) 

W = D 2 W = D 3 W= 0 = 0 at z = 0 and 1. 

(24) 

The case of two free boundaries, 
though a little artificial is the most 
appropriate for stellar atmospheres. Using 
the boundary conditions (24), we can show 
that all the even order derivatives of W must 
vanish for z = 0 and z = 1 and hence the 
proper solution of W characterizing the 
lowest mode is 

W = W 0 sin nz ; W 0 is a constant. (25) 

Substituting equation (24) in (22), we obtain 
the dispersion relation 



R^xP = 



1 + 



— n 2 F io- 



(1 + x) 



(1 + x + EfriaJ + D Al (l + x) 2 
(1 + x + E^ia-J, 



(26) 



x = 



where R 1 = - , D Al = - A , 
io t = ^,P = n 2 P t . 

Equation (26) is required dispersion relation 
accounting for the onset of thermal 
convection in Walters' (Model B ) elastico- 
viscous fluid in a Brinkman porous medium. 

4. STABILITY OF THE SYSTEM 
AND OSCILLATORY MODES 

Here we examine the possibility of 
oscillatory modes, if any, in Walters (Model 
B') elastico-viscous fluid due to the presence 
of viscoelasticity, medium perpeability and 
gravity field. Multiply equation (21) by W 
the complex conjugate of W, integrating 
over the range of z and making use of 
equations (22) with the help of boundary 
conditions (24), we obtain 



+ 



[ 1 + g_ F ) ff ] 7l _^ (/a 
E Pl G*I 3 )-D A (I 4 + a 2 I 5 -) = 0, (27) 
where 

/i = fi(\DW\ 2 + a 2 \W\ 2 )dz, 

1 2 = C(\D8\ 2 + a 2 \e\ 2 )dz, 

1 3 = ti\0\ 2 dz, 
U = tf\DW\ 4 dz, 

h = ti(2\DW\ 2 + a 2 \W\ 2 )dz. 

The integral part Ij-I 5 are all positive 
definite. Putting a = ia t in equation (26), 
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where a t is real and equating the imaginary 
parts, we obtain 



= 0, (28) 



Equation (28) implies that a t = 0 or 
Oi ^ 0 which mean that modes may be non 
oscillatory or oscillatory. The oscillatory 
modes introduced due to presence of 
viscoelasticity, gravity field and medium 
permeability. 

5. THE STATIONARY 
CONVECTION 

For stationary convection putting er = 0 in 
equation (25) reduces it to 
(l+x) 2 



Ri = 



xP 



•[l + (l + x)D Al ] 



(29) 



which expresses the modified Rayleigh 
number R[ as a function of the 
dimensionless wave number x and the 
parameters D A , P and Walters' (Model B ) 
elastico-viscous fluid behave like an 



ordinary Newtonian fluid since elastico- 
viscous parameter F vanishes with a . 

To study the effects of Darcy 

number and medium permeability, we 

dRi , di?! 



examine the behavior of 
analytically. 

From equation (29), we get 



and 



dP 



dR t 
dD & 



(1+x) 3 
xP 



(30) 



which shows that Darcy number 
stabilizing effect on the system. 
It is evident from equation (29) that 
dR t (1+x) 2 
dP ~ xP 2 



[l + (l + x)D A J 



has 



(31) 



From equation (30), we observe that 
medium permeability has destabilizing 
effect. 

The dispersion relation (29) is 
analyzed numerically. Graphs have been 
plotted by giving some numerical values to 
the parameters, to depict the stability 
characteristics. 




Fig.l. Variation of Rayleigh number R 1 with Darcy number D Al for P = 2 for fixed wave numbers x = 0.2, 
x = 0.5 and x = 0.8. 
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Fig.2. Variation of Rayleigh number R 1 with medium permeability P for D Al 
x = 0.2, x = 0.5 and x = 0.8. 



5 for fixed wave numbers 



In fig.l, Rayleigh number ft x is 
plotted against rotation D A for B = 3, P = 2 
for fixed wave numbers x = 0.2, x = 0.5 and 
x = 0.8. This shows that Darcy number has a 
stabilizing effect for fixed wave numbers x = 
0.2, x = 0.5 and x = 0.8. 

In fig.2, Rayleigh number R t is 
plotted against medium permeability P for 
D A = 5 for fixed wave numbers x = 0.2, x 
= 0.5 and x = 0.8. This shows that medium 
permeability has a destabilizing effect. 

6. CONCLUSION 

The onset of thermal convection in 
Walters' (Model B ) elastico-viscous fluid 
in a Darcy-Brinkman porous medium has 
been investigated. The dispersion relation, 
including the effects of Darcy number, 
medium permeability and viscoelasticity on 
the thermal instability in Walters' (Model 
B ) fluid is derived. From the analysis, the 
main conclusions are as follows: 



(i) For the case of stationary convection, 
Rivlin-Ericksen elastico-viscous fluid 
behaves like an ordinary Newtonian 
fluid as elastico-viscous parameter F 
vanishes with a . 

(ii) The expressions for dRl and are 

v ' v dD Al dP 

examined analytically and it has been 
found that the Darcy number has 
stabilizing effect whereas the medium 
permeability has a destabilizing effect 
on the system. The effects of Darcy 
number and medium permeability on 
thermal convection have also been 
shown graphically in figures 1 and 2. 

(iii) The oscillatory modes introduced due to 
presence of viscoelasticity, gravity field, 
and medium permeability. 
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